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Abstract. 

We briefly review the SUGRA, GMSB and AMSB supersymmetry breaking models. 
We then discuss the phenomenological differences between them and consequent 
characteristic experimental signatures. This is followed by a review of the discovery 
potential for supersymmetry at the Tevatron, LHC and a future e + e _ linear collider. 



1. Introduction 

In the past twenty years there has been a great deal of theoretical study of low energy 
supersymmetry and a number of experimental searches. However, these searches have 
found no evidence for supersymmetry. Within the next ten years a number of new 
collider experiments will probe higher energies and low energy supersymmetry will either 
be discovered experimentally or will no longer be relevant to the problem of electroweak 
symmetry breaking. 

After briefly surveying the models which are used in experimental studies we will 
discuss the discovery potential and signatures for current and future experiments. 



2. SUSY Models 



We will consider the Minimal Supersymmetric Standard Model (MSSM). This is the 
supersymmetric extension of the the Standard Model (SM) which has minimal particle 
content, i.e. it contains two Higgs doublets rather than the one of the SM and the 
superpartners of the Standard Model fields. This model has one less parameter than 
the Standard Model, provided that supersymmetry (SUSY) is unbroken. There is an 
additional parameter /i which gives mixing between the two Higgs doublets but the 
couplings in the scalar potential are constrained by SUSY. However, as the superpartners 
have not been observed, SUSY must be broken in such a way as to not reintroduce the 
quadratic dependence of the Higgs mass on the cutoff scale. This is called soft SUSY 
breaking and leads to a large number of additional parameters: soft SUSY breaking 
masses for the scalars and gauginos; soft SUSY breaking A terms, which couple two 
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sfermions and a Higgs; and B terms, which couple two Higgs bosons. The trilinear 
A terms are only important for the third generation sfermions as they enter in terms 
proportional to the fermion masses. The Z boson mass is given in terms of the other 
parameters of the model. 

In the MSSM R-parity related to baryon (B) and lepton number (L) by 
Rp = (_i) 3B + 2S, + L ; i s assumed to be conserved which means that SUSY particles can 
only be produced in pairs and the lightest SUSY particle (LSP) is stable. As the LSP is 
stable it must be neutral as charged stable matter is excluded by cosmological arguments. 
This implies that all SUSY events have two LSPs in them leading to missing transverse 
energy in SUSY events because the LSPs do not interact in the detector. R-parity need 
not be conserved and this leads to very different experimental signatures jl]. 

It is difficult to use the MSSM for detailed experimental studies due to the large 
number of free parameters it contains. Furthermore large ranges of the parameters are 
excluded and a real model will surely have far fewer fundamental parameters. Therefore 
models of SUSY breaking are used which are motivated by theoretical ideas. All have 
a common feature; SUSY is broken in some hidden sector and then transmitted to the 
MSSM fields. The models differ in how this transmission is accomplished: 

SUGRA In supergravity models || all the scalar masses (M ), the gaugino masses 
(Mi/ 2 ), the A and B parameters are assumed to be unified at the grand unified 
(GUT) scale (~ 10 15 GeV). As the model predicts Mz in terms of the other 
parameters it is possible to use tan/3 = t>i/t>2, the ratio of the vacuum expectation 
values (VEVs) for the two Higgs doublets, and the known value of Mi to fix B 
and | /x | . This leaves five parameters Mo, M1/2, A, sgn/i, tan/3 which completely 
determine the mass spectrum and decay patterns of the particles. The gluino 
mass (M|) is strongly correlated with M1/2 and the slepton mass with M . The 
LHC experiments have defined several SUGRA points which are often used in 
simulations ||[4|]. 

GMSB The Gauge Mediated SUSY Breaking model |J aims to solve the flavour 
changing neutral current problem which is generically present in SUGRA models by 
using gauge interactions rather than gravity to transmit the SUSY breaking. The 
messenger sector consists of some particles, X, which have SM interactions and are 
aware of SUSY breaking. The simplest choice is to have the messenger particles 
in complete SU(5) 5 or 10 representations to preserve the GUT symmetry. The 
fundamental SUSY breaking scale F must be such that yF < 10 10 GeV or SUGRA 
breaking will dominate. The gaugino masses occur at one-loop while the squark 
and slepton masses occur at two-loop. The LSP is an almost massless gravitino so 
that the sparticles decay as in a SUGRA model followed by the decay of the next- 
to-lightest SUSY particle (NLSP) to the gravitino. The NLSP need not be neutral 
and its lifetime is model dependent. This model has six parameters: A = F/M the 
scale for the SUSY masses; M > A the messenger mass scale; N 5 > 1 the number 
of 5 + 5 messenger fields; the ratio of the vacuum expectation values tan/3; sgn/x 
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Figure 1. Examples of mass spectra in SUGRA, GMSB and AMSB models for 
. The other parameters are: M = 100 GeV, M 1/2 = 200 GeV 



tan /3 = 3 and sgn // = 

for the SUGRA model; M = 100 TcV, N 5 = 1, A = 70TeV for the GMSB model; 
Mo = 200 GeV, M 3/2 = 35 TcV for the AMSB model. Plot taken from §. 



the sign of the n parameter; and C grav > 1 which controls the NLSP lifetime. 

AMSB The super-conformal anomaly || is always present and predicts sparticle masses 
in terms of M 3 / 2 , the gravitino mass. The simplest version of this model predicts 
tachyonic sleptons and therefore some other SUSY breaking mechanism must be 
present in order to get a realistic spectrum. One way to do this is to add a universal 
scalar mass (mAMSB) or new very heavy fields (DAMSB). The (mAMSB) model 
has four parameters: M the universal scalar mass; M 3 / 2 the gravitino mass; tan/3 
the ratio of Higgs VEVs; and sgn fi. The DAMSB model has five parameters: M 
the mass of the new fields; n the number of new fields; M3/2; tan/3; and sgn/x. 
The AMSB model has one important feature in that the LSP is mainly wino like 
and almost degenerate with the lightest chargino, xt 1 but this feature is lost in 
DAMSB. 

All of these models are implemented in the ISA JET event generator [[?]]. 

The spectra of the SUSY particles in these models can be very similar, an example 
of each is shown in Fig. [1]. The main differences in these models are the ratios of the 
squark to slepton masses and the differences between the electroweak gaugino and gluino 
masses. In the SUGRA model the scalar masses are universal at the GUT scale and 
therefore the differences in the masses come from the renormalization group evolution 
(RGE) to the electroweak scale. However, in the GMSB models the masses at the SUSY 
breaking scale are proportional to the relevant gauge couplings and therefore the strongly 
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interacting squarks are heavier than the sleptons, relative to the SUGRA model. In the 
AMSB model the sleptons are very light compared to the squarks unless the universal 
scalar mass is very large. Similarly for the gauginos the splitting of the gluino and 
electroweak gaugino masses in the SUGRA model comes from the RGE however in the 
GMSB model the gaugino masses are proportional to the gauge couplings at the SUSY 
breaking scale and therefore the strongly interacting gluino is heavier than the weakly 
interacting gauginos. In the AMSB model the soft breaking mass for the gluino is larger 
than the soft breaking masses for the electroweak gauginos at the GUT scale which 
gives a bigger splitting between the gluino and electroweak gauginos masses than in 
SUGRA models. The nature of the lightest neutralino is also different in the different 
models. In the SUGRA model it is usually mainly bino whereas in the AMSB model 
it is mainly wino and degenerate with the lightest chargino. As discussed above, the 
lightest neutralino in the GMSB model is the Gravitino and the NLSP, which can be 
neutral or charged, is most important for phenomenology. 

Obviously once supersymmetry is discovered it will be important to make accurate 
measurements of SUSY particles masses and couplings in order to investigate the 
model of SUSY breaking. However at present, given we have seen no evidence for 
supersymmetry in experimental studies, the main interest is in simulating models which 
give qualitatively different experimental signatures so that we can be certain that all 
variants can be observed. For example the DAMSB model is very similar to a SUGRA 
model and has not therefore been subjected to many detailed studies, whereas in the 
mAMSB model which has an almost degenerate lightest neutralino and chargino the 
dominant chargino decay mode is xl ~ * 7l+ Xi which is very different from SUGRA 
models and therefore of more interest. 

There are a number of signatures which are characteristic of supersymmetry, 
regardless of the model of SUSY breaking: missing transverse energy; a high multiplicity 
of high transverse momentum jets; many isolated leptons; copious b-jet production; a 
large rate of Higgs production; isolated photons; and quasi-stable charged particles. It 
should be noted that not all of these signals are present in all models and that production 
of any heavy object will give some of these signals. 

In order to simulate supersymmetry it is essential to have a consistent model. 
We cannot consider one sparticle in isolation because all the supersymmetric particles 
which are kinematically allowed will be produced. In hadron colliders production of the 
squarks and gluinos dominates provided the centre-of-mass energy is high enough. The 
production of those sparticles which only have electro- weak couplings may be dominated 
by the decays of squarks and not by direct production. The dominant backgrounds 
at the LHC are combinatorial from SUSY events after some simple cuts are applied. 
The situation at the Tevatron where electroweak gaugino production may dominate 
and Standard Model processes are the most important source of background is very 
different. At a lepton collider where the full spectrum is unlikely to be accessible and 
the beam energy and polarization can be used to separate sparticles the situation is 
much simpler. 




Figure 2. Tevatron reach in the Tri-lepton channel in the Mq,M\iq, plane, for fixed 
values of Aq = 0, /j > and (a) tan/3 = 5 or (b) tan/3 = 35. Results are shown for 2, 
10 and 30 fb _1 integrated luminosity. The curves require the observation of at least 5 
events and are 3<7 exclusion contours. The cross-hatched region is excluded by current 
limits on the superpartner masses and the dot-dashed lines correspond to the projected 
LEP-II reach for the chargino and lightest Higgs masses. Figure taken from |9j,|To[] . 



3. Tevatron 

In Run I of the Tevatron neither CDF or DO claimed the discovery of any signal for 



super symmetry A number of limits were obtained which we will not discuss here [|1"2|,|T3 
Run II of the Tevatron which has both an increase in the centre-of-mass energy to 2 TeV 
and one hundred times the luminosity will extend the mass range but the search reach 
is limited. Due to the centre-of-mass energy the production of squarks and gluinos 
may not dominate and the best channel for the discovery of SUSY may be gaugino 
production, i.e. production of X2X1 ^ + ^~X?^ +Z/ Xi- Fig-H shows the search potential 
in this channel at Run II of the Tevatron for different integrated luminosities. The 
background to this process is dominated by gauge boson pair production followed by 
leptonic decays with WZ/7 being the dominant background process. 

If this signal is observed structure in the £ + £~ mass distribution will constrain the 
X\ and X2 masses, this will be discussed in more detail for the LHC in the next section. 
However as can be seen from Fig. |3] for those points which are still allowed by the LEP 
limits the endpoint of the distribution which gives information on the mass difference 
of the lightest two neutralinos will be difficult to measure. 

It is possible to extend the search reach by using channels involving jets and missing 
transverse energy. If supersymmetry is discovered at the Tevatron it will determine the 
mass scale of some of the particles. The LHC will then make detailed measurements of 
the SUSY spectrum. 
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Figure 3. Opposite sign, same flavour dilepton mass reconstruction for five study 
points given below and the WZ and tt background. The background is included in 
all the histograms. It should be noted that since this plot was produced point 1 has 
been ruled out by searches at LEP. Point 1 has M = 100 GeV, M 1/2 = 200 GeV, 
A = OGeV, tan/3 = 3, sgn/i = +; point 2 has M = 140 GeV, M 1/2 = 175 GeV, 
A = OGeV, tan/3 = 35, sgn^ = +; point 3 has M = 200 GeV, M 1/2 = 140 GeV, 
A = -500 GeV, tan/3 = 35, sgn/i = +; point 4 has M = 250 GeV, M 1/2 = 150 GeV, 
A = -600 GeV, tan/3 = 3, sgn/i = +; point 5 has M = 150 GeV, M 1/2 = 300 GeV, 
Ao = OGeV, tan/3 = 30, sgn/i = — and non- universal GUT-scale Higgs masses 
M Hu h 2 = 500 GeV. Plot taken from |Q. 
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Figure 4. Discovery potential of the LHC for final states with at least two jets, 
> 100 GeV and at least one isolated lepton. An integrated luminosity of 100 fb _ 
has been assumed. Plot taken from JLJ] . 



4. LHC 

Many detailed studies of the search potential of the LHC have been performed by both 
the ATLAS || and CMS Q collaborations. At the LHC the strongly interacting squarks 
and gluinos are predominantly produced unless they are very heavy giving large SUSY 
production cross sections. This allows hard cuts to be applied to the data in order 
to reduce the Standard Model background leaving combinatorial backgrounds from the 
SUSY events themselves as the dominant source of background. Most of the studies 
have assumed the full integrated luminosity of the LHC which enables hard cuts to be 
used. In practice the discovery of SUSY will probably take less time and use weaker 
cuts in which case the Standard Model background will be more important. 

The LHC studies can be broadly grouped into two categories: the first attempts 
to find some signal indicative of SUSY and find an excess in order to discover 
supersymmetry by exploiting inclusive signatures; the second then tries to make use 
of more information from the SUSY events in order to measures masses and other 
parameters of the model. 

There is a very large range of accessible masses in inclusive signals, i.e. jets, 
leptons and missing transverse energy (Jftr). Fig.|] shows the mass reach in SUGRA 



models for the CMS experiment [ 14| . This covers all the interesting theoretical range, 
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Figure 5. (a) The reconstructed bb mass distribution for events passing h — > bb 
selection cuts. The distributions are for the SM background (shaded), the total 
SUSY+SM background (dashed) and the summed signal+background for SUGRA 
point 1. (b) The 5cr discovery contour curves for h — > bb from SUSY cascade 
decays in the Mq,M 1 / 2 plane for tan/3 = 10, sgn/i = +. The expected number of 
reconstructed h — > bb events are also shown. The dark shaded regions are theoretical 
or experimentally excluded. Both (a) and (b) assume 30 fb -1 integrated luminosity 
and are taken from B . 



i.e. M| < 2.5 TeV. It is useful to define global variables for SUSY searches. For example 
if we consider events with at least four jets and missing transverse energy the variable 

M eS = Pt,i + PT,2 + PT,3 + Pta + $t, (1) 

where pr,i is the transverse momentum of the ith jet, is very useful. The peak in the 
M cS distribution correlates well with the SUSY mass scale where M SU sy = min(M fi , M g ). 



This can determine the squark/gluino masses to about 15% |p"5* . 

There have been many studies of techniques to reconstruct sparticle masses and 
properties. Here we will illustrate the techniques by choosing examples from case 
studies. In general both the squarks and gluino are produced. Then depending on the 
relative masses one decays into the other as these decays occur via the strong interaction. 
The weak gauginos are then produced in the decays of the lighter strongly interacting 
particles, for example q^ — > x°q. In most models a significant number of second-to- 
lightest neutralinos are produced. These neutralinos then either decay via x\ ~* Xi n 
or x\ — > Xi^ + £~ i possibly via either an intermediate slepton \\ £ + £~ — > x\t r t~ 
or Z boson \\ ~* ^Xi ~^ Xi^ + ^~ ■ The Higgs decay mode tends to dominate if it is 
kinematically accessible. Most studies have used these decays as a starting point for 
mass measurements. Many other SUSY particles can then be identified by adding more 
jets or leptons to reconstruct other particles in the decay chain. 

If the decay of the x\ — > x?h exists then approximately 20 % of SUSY events contain 
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Figure 6. 



(a) Dilepton X ° ~> ~* signal at SUGRA point 5 

for 100 fb -1 including backgrounds after flavour subtraction, i.e. the plot shows 
e+e~ + - e ± fi T . (b) Dilepton x° — * x\^^ distribution for SUGRA point 

4 (solid) and SM background (shaded) with 30 fb -1 . Both plots are taken from H]. 



h — > bb. In these models the Higgs would be discovered in SUSY events at the LHC 
rather than by the traditional Standard Model-like Higgs searches. The mass of bb 
pairs is shown in Fig. 0(a) after the following cuts: $t > 300 GeV; more than two jets 
with pt > 100 GeV and more than one with \rj\ < 2; no isolated leptons to suppress 
the tt background; only two b-jets with p T h > 55 GeV and \rj\ < 2; AR h ^ < 1.0 again 
to suppress the tt background. This gives a clear peak in the bb mass distribution at 
the Higgs mass. The SM background is very small and the dominant background is 
from other SUSY decays. This method works over a large region of parameter space in 
SUGRA Models, Fig.|(b). 

In the regions of parameter space where the decay mode — > x?h is n °t 
kinematically accessible, the reconstruction of the leptonic decay mode is necessary 
in order to constrain the m ass. The important decay modes are either via a real 
slepton, Fig. ^(a), or via virtual sleptons and gauge bosons Fig. §|(b). In Fig.|6|b there is 
also a peak at the Z mass due to the production of Z bosons in other SUSY decays. 

The leptons produced in neutralino decays can be combined with the other decay 
products of the heavier SUSY particles in order to reconstruct them. For example the 
X2 is often produced in left squark decays, qx — ► qx° - * q^ + ^~Xi- m order 

to identify this decay chain the following cuts are applied: two isolated leptons with 
transverse momentum, p? > 10 GeV; more than four jets one with p? > 100 GeV and 
the rest with p? > 50 GeV; missing transverse energy Ifi? > max(100.0, 0.2M e ff). The 
mass of the q££ system has a maximum at 
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Figure 7. (a) Distribution of the larger t + i~q mass for M u > M$ ax /\/2. (b) The 
smaller of the two £ + £~q masses for the signal, (c) The £ ± q mass distribution for 
combinations with M^+^ q < 600 GeV. Plots taken from ||. 



where M? L is the squark mass and M^o 2 are the lightest and next-to-lightest neutralino 
masses respectively. This distribution is shown in Fig.[7|(a). The minimum of this 
distribution, Fig. 0(b), and the £q distribution, Fig. 0(c), also provide useful information. 
This system has four constraints, i.e. the upper edges in £q, ££ and ££q distributions, 
and the lower edge in the ££q distribution, and four unknowns, i.e. M^o, M^o, M? L 
and M~ which can therefore be reconstructed. The errors on the x®, x%i and £r 
masses are 12%, 6%, 3% and 9%, respectively. The mass of the unobserved LSP is 
determined, Fig.|8]. The errors on the particle masses are strongly correlated and a 
precise determination of one mass would reduce the errors on the rest [16| . 

In GMSB models a similar type of analysis is possible. For example, consider a 
case where the NLSP is the lightest neutralino which then decays to a gravitino and a 
photon, i.e. x\ — * 7G. As in SUGRA models the x° can decay leptonically to give the 
lightest neutralino. This gives the decay chain, x® ~~ * ^ + ^~Xi ~^ ^ + ^~7G. If we require 
M eS > 400 GeV; $t > 0.1M cff ; at least two leptons and two photons where photons 
and electrons have px > 20 GeV and muons px > 5 GeV. Here information can be 
obtained from the £ + £~, £ + £~^i and £ ±r y mass distributions. There are two structures 
in the £ ±r ) distribution and therefore these distributions provide four constraints which 
are sufficient to measure the xli X2 an d £r masses in this model. 

In AMSB models the signatures are very similar to SUGRA models except that the 
Xi may not be observable due to the small mass difference between the lightest chargino 
and the LSP. 

The signatures of SUSY in models with R-parity violation are very different because 
the LSP can decay inside the detector which means there may be no missing transverse 
energy. The LSP will decay to give either leptons ||[17|], leptons and jets or just 
jets P,|18[. In these models the LSP mass can be found by reconstructing all its decay 
products which often allows the LSP mass to be measured with greater precision than 
is possible in SUGRA models. 
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Figure 8. (a) Distribution of the fractional difference between the reconstructed 
and true squark mass, (b) Distribution of the fractional difference between the 
reconstructed and true %5 mass. Both plots were taken from ||. 



5. Lepton Colliders 

Any future linear collider will start operation after the LHC and therefore the study 
of supersymmetry at such a machine will be in the context of what we already know 
from the LHC. A collider with centre-of-mass energy of less than 1 TeV will probably 
concentrate on the sparticles which only interact weakly, i.e. the electroweak gauginos 
and sleptons. A higher energy collider may also be able to produce the squarks, however 
at any lepton collider production of gluinos is difficult, unless they are lighter than the 
squarks in which case they are produced in squark decays. 

An e + e~ machine has a number of advantages over a hadron collider: the energy 
of the beam provides a kinematic constraint; polarization of the electrons and possibly 
positrons provides a powerful tool and the signal to background ratio is of order one 
before cuts are applied. The great power of such a facility is the ability to make precise 
measurements of the masses and couplings of sparticles which may already have been 
observed at the LHC. All the examples we will consider are taken from M. Similar 
studies can also be found in [[19] . 



If we consider, for example, the production of smuons, e + e~ — > /a + /C — > ~> 
the events will have a pair and missing energy. The energy spectrum of the muons 

is shown Fig.|9|. This spectrum is flat apart from beamstrahlung, initial-state radiation 
and resolution effects at the high edge. The end points of this distribution can be related 
to the masses of the smuon and lightest neutralino. Using this process both the x? and 
jl masses can be determined to ~ 0.5 %. If the polarization of the beams is changed 
the amount of left and right slepton in the mass eigenstate can also be determined. In 
the case shown in Fig. ||] the machine is below threshold for the other SUSY particles, 
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Figure 9. The energy spectrum of the muons produced in the process 
e + e~ — > /}r/}r — > /^ _ XiM + Xi & t \/s = 320 GeV for an integrated luminosity of 160 fb _1 . 
Plot taken from ||. 



apart from x^x^ an d therefore the SUSY background is small. The Standard Model 
background is even smaller after event selection. 

The heavier gauginos can also be produced, e.g. e + e~ — ► X2X2 ~~ ^ ^ + ^~^ + ^~XiXi 
at sufficiently high centre-of-mass energy. The dilepton masses and energies for this 
process are shown in Fig.[l(| The masses of both the Xi an d X2 can be determined from 
the di-lepton energy spectrum with typical errors of about 0.2%. The di-lepton mass 
spectrum gives additional information on the mass difference between the two leptons 
and this difference can typical be measured with a precision of better than 50 MeV. 

These mass measurements will be more accurate than those achievable by the LHC, 
if the lepton collider has sufficient energy to produced the particles. Furthermore it 
should be possible to determine the mixings in the gaugino sector which is difficult at 
the LHC. Note that, as discussed above, the extraction of masses at the LHC may result 
in strongly correlated errors. Measurements from a linear collider of some masses could 
therefore be used to improve the errors on the LHC results for the heavier particles, e.g. 
the squarks, which the lepton collider may not be able to produce. 

6. Conclusions 

The upgraded Tevatron has some search potential for supersymmetry however, given the 
current limits, the event rates will be low. The most promising signal is from tri-leptons 
produced following the production of xf X2, but the search range is limited. 

The production rate of SUSY particles at the LHC is very large unless the squarks 
and gluinos are very heavy. Looking for signals of direct production of gauginos and 
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Figure 10. (a) The di-lepton mass and (b) di-lepton energy spectrum for 
e + e~ -> X2X2 e+i-j^e+e-xl at = 320 GeV for an integrated luminosity of 
160 fb" 1 . Plot taken from §|. 

sleptons is difficult but not impossible (jet vetoes have to be used). It is very difficult to 
observe the heavier gauginos unless they are strongly mixed because otherwise they are 
mainly Higgsino and therefore do not couple to the squarks and so cannot be produced 
in their decay. 

A future e + e~ collider would be a very powerful tool for precise measurements 
of masses and couplings provided that the energy is high enough. A few precise 
measurements made with such a machine could, in combination with the LHC 
measurements, greatly constrain the underlying model. 

We are approaching the end of an era. Low energy supersymmetry has been studied 
for the last twenty years, without any experimental verification. However the within 
the next eight years or so, the searches will reach high enough mass scales so that it will 
either be discovered or cease to be relevant to the problem of electro-weak symmetry 
breaking. 
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